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Abstract. Semi-analytical model for calculating acoustic response to a line-focused 
laser pulse in an optically absorptive isotropic cylinder is proposed and implemented. 
It is assumed that the laser input is absorbed over the volume and thus creates a radially 
distributed thermoelastic source. Closed-form solution is obtained in the Fourier domain. 
Two inverse Fourier transforms in frequency and circumferential wave number yield the 
sought waveforms of acoustic response in the time-space domain. Numerical simulation 
is compared to the calculation based on a surface dipole source. The two signals have 
essentially different shapes of the wave-arrival peaks when the optical penetration is large 
enough. 
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1. Introduction 

The laser ultrasonics technique is a unique tool for acoustic- wave generation in objects with 
a curved surface. Much work has been done by this means for generating surface and bulk 
acoustic waves in opaque cylindrical targets with a view of their non-destructive testing 
and evaluation, see recent publications [1-4] and the bibliography therein. There exists a 
well developed theory for the case when, at a given laser wavelength, the target (either 
cylindrical or planar) is opaque and hence the source is located at the surface. Given so, 
the ablation model applies for high enough local energy such that causes ablation at the 
laser-heated region of the surface, and the dipole (Scruby's) model applies for relatively 
low incoming energy [5]. In this latter, thermoelastic regime, the dilatation of an infinitely 



small volume adjacent to the surface produces a surface source equivalent to a force dipole. 
The dipole strength has been derived in [6, 7]. 

Thermoelastic regime for the optically absorptive materials implies that the laser input 
penetrates into a target and creates a volume distribution of sources therein. Hence the 
assumption of a surface source is no longer appropriate. Laser generation of acoustic waves 
in optically absorptive materials has been considered for the case of a planar target (plate) 
[8]. The next natural step is extending laser ultrasonics treatment to optically absorptive 
cylindrical targets. To the best of the authors' knowledge, such a case has not yet been 
studied. The present paper proposes a two-dimensional (2D) semi-analytical model for 
calculating acoustic response to a linc-focuscd laser pulse in an optically absorptive isotropic 
cylinder. Thermoelastic regime of generation is assumed and a typical case of materials with 
negligible thermal conductivity is treated. Closed-form solution is obtained in the Fourier 
domain to the inhomogeneous ordinary differential system with a radially distributed 
thermoelastic source and with the traction-free boundary condition at the cylinder surface. 
Performing two inverse Fourier transforms yields the sought waveforms of acoustic response 
in the time-space domain. This calculation is carried out in the paper for two different 
source settings. One of them formally treats the case of an opaque cylinder by assuming a 
vanishingly small optical penetration or, which is equivalent, a buried source tending to the 
cylinder surface. The other considers the actual case of an optically absorptive cylinder, 
where the exponentially decaying law of optical absorption produces a distributed source 
with large enough radial penetration length. The results obtained for these two settings 
are compared with the calculation based on Scruby's model of a surface dipole source. It 
is observed that the waveforms, computed for a buried source tending to the surface, agree 
very well with those following for Scruby's surface-source model, whereas the waveforms, 
computed for a distributed source, have essentially different shapes comparatively to the 
prediction following for the surface source. The former observation serves as a test for our 
model, while the latter highlights the effect of laser-energy penetration into the bulk of an 
optically absorptive cylinder. 

2. Background 

Consider an infinite homogenous isotropic elastic cylinder of radius R, whose surface is 
heated by the line-focused laser pulse independent of the axial coordinate z. The laser 
energy is supposed to be low enough to preclude ablation. It is assumed that viscoelasticity, 
thermal conduction and the contribution of deformation to the entropy change may be 
disregarded. The thermoelastic generation of acoustic wave is therefore described by the 
following simplified system of governing equations: 

V • (7 = pii; = Atr (e) I + 2 ^e-/3TI; pcpf = q, (1) 

where a is the stress tensor including mechanical and thermal stresses, e is the strain tensor, 
u is the displacement, p is the density, A and p are the Lame coefficients, /3 = (3A -|- 2;u) ar 
is the thermal modulus (ar is the dilatation coefficient), T is the temperature rise, Cp is 
the specific heat, q is the power density of the line laser source, I is the identity matrix 
and tr denotes the trace. 



The heat input q{r,6,t) creates a volume distribution of thermal sources inside the 
cylinder through the optical absorption. Assuming normal incidence of the Gaussian laser 
beam [9] and the exponential fall off of the incoming energy E, it follows that 

q{r,e,t)= aEgie)6it)e-''^''-^\ (2) 

where 

1 /41n2 -(4in2)4 o . ^ 
g(e) = -\ 7 = 2 arctan— (3) 
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describes the angular extent of the beam with width b. By (1)3 and (2), 

T(r,0,t) = — 5Wff(t)e-"(«-'-), (4) 

pCp 

where H is the Heaviside function. The temperature field (4) specifies the thermal stress 
C7(^) = -pT in (1). 

Thus, the thermoelastic problem is reduced to solving the 2D plane-strain equation 
of motion for the displacement u(r, 9, t) = UrBr + ugeg in the cylinder cross-section 
{er,eg). Divergence of the thermal stress plays the role of a source term. The equation 
is complemented by the traction-free boundary condition (BC) e^cr = at the cylinder 
surface r = R. 



3. Solution in the Fourier domain 

3.1. Setup 

Consider the above 2D problem in the Fourier domain {r, i^, w}, where v is the 
circumferential wave number and uj the angular frequency. Thereafter in this section the 
factor e^i''^-'^^) is omitted, the same notations for the transformed functions are used, and 
their dependence on 1/,^ is understood and suppressed. According to (1) and (4), we write 

a (r) = (r) + a^^) (r) I, 
^(m) ^ xtv [e (r)] I + 2 (r) , cj(^) (r) = -/JToe-^^-^"'') , ^ ' 

where Tq a) involves double transform of (4) in 6 and t. 

The equation of motion in the Fourier domain reduces to the ordinary differential system 
in u (r) with the source term provided by Vcr^-^^ (r) . Using Helmholtz decomposition, this 
equation may be re-cast in the form of two uncoupled Bessel equations for the longitudinal 
and shear elastic potentials (p (r) and ip (r), 
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where k\ = pu? j (A + , k\ = pu? j p and h = (3Tq/ (A + 2p) . The traction-free BC, 
which couples (p (r) and V' {t) , is 
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For the future use, introduce the fohowing abridged notations of the values taken at the 
boundary r = R: 



Kl,t = kL,TR, Bt,l = 
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dx 



(9) 



Note that expressing the results in terms of Bt,l is helpful for numerical implementation, 
because the exponential growth of the Bessel functions of the order u S> Kl^t is avoided 
this way (see [10]). 



3.2. Green function 

First, we solve an auxiliary problem for a buried point source at ro (0 < ro < R) under the 
homogeneous BC on the mechanical stresses (see (8)i), 
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= 0, 
(10) 

where the factor R^ in (10)i is added to merely keep physical dimensions intact. The 
partial solution to the given inhomogeneous Bessel equation on $ (r|ro) is [11] 

g (rjro) = bR^^ [J, {kLv) (kLVo) H (ro - r) + {kLv) J. (fc^ro) H {r - ro)] , (11) 

where and Yn are the Bessel functions of the first and second kind. Thus the solution 
to (10) is sought in the form 



$ (r|ro) 
* (rjro) 



g (rjro) + C$ (ro) J„ (kLv) 
(ro) J„ (krr) 



(12) 



where Jjy (kir) and Jjy {kxr) are the solutions of the homogeneous Bessel equation on the 
potentials. Substituting (12) into the BC of (10) and appealing to (8)2 provides an algebraic 
system of equations on the coefficients C$ (ro) and Cj. (ro) . It has the form 
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with 
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where the use of (9) has been made, and B^^^ differs from hy replacing with Yi^. 
Note that the matrix (first-rank tensor) A results in (13) from applying the boundary stress 
operator S (i?) to the solutions of the homogeneous equation. Hence D = 0, where 

is the dispersion equation for the eigenmodes of the traction-free cylinder. By (13), 

Inserting (16) along with (11) into (12) completes finding the solution $ (r\ro) and ^' (r|ro) 
to (10). 

3.3. Displacement for a distributed source 

Once the solution $(r|ro) and ^ (r|ro) to (10) has been found, it remains to take its 
convolution with the radially distributed source ~ Q-a{R-r) Jqj. obtaining the solution 
(f (r) and V' (f) to Eq. (6), and to invoke Hclmholtz decomposition u = V(/3-|- V x (tpe^) for 
finding u (r) . Because these two are linear operations, they may be applied in the optional 
order. Thus 

u (r) = £ {V$ (r|ro) + Vx[^ (r|ro) e,]} e-"(^-'-°Vodro, (17) 

where G(r|ro) = V$(r|ro) + V x [^'(r|ro)e^] is the displacement response to a buried 
source Vcr^-^^ (r) ~ yi(l^_Llll^ Note that application of the boundary-stress differential 
operator to g {r\rQ) (see (12)) in the integrand of (17) taken with r, tq R provides 
an extra delta-function term, which compensates (T^^-* (R) = —^Tq in the BC (7) for the 
distributed source. 

In the context of laser ultrasonics, the important particular case is the displacement at 
the surface r = i2. It can be written in the following closed form: 

u (B) = f ■) = ( -f], (18) 

where 
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1= C J, {Klx) e-"^(i-^)xdx {x = r/R) , C/q = (19) 
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and D is given in (15). 

Numerical simulations in §4.2 involve the model case of a buried source located at 
ro = R — e tending to the surface. For this case, u (R) = -^G {R\R — e) is given by (18) 
with I = Jy [Kl) (the equality is to within arbitrarily small e/R). Obviously, this case 
may also be seen as a limit of vanishingly small optical penetration length a~^. 



4. Numerical treatment and the simulation results 

4-1- Inverse Fourier transform 

2D Fourier transform, discrete in circumferential number v and continuous in frequency to, 
has next been performed to obtain the sought acoustic response in the time-space domain 
{r, 9, t}. The calculation has been made for the radial component of response at the surface, 

which is most usually the quantity measured by the laser ultrasonics techniqiic. Denote its 
abovc-dcrivcd value in the transform domain as Ur (R) = Ur {R, ^) ■, where {R) is given 
by Eq.(18). Fast Fourier Transform algorithm in frequency has been applied, with a small 
imaginary part added to w {uj* = to—jS where S > 0, e.g. [12]) in order to avoid zeros of the 
dispersion equation D = 0, see (19) and the remark to (15). On taking Fourier transform 
in v, the use has been made of the fact that, by (18), tir(r, z^, w) is an even function of v. 
As a result the sought response Ur{R, 9,t) is described by the formula 

/-oo C ~ ^ 

Ur{R, e, t) = / <^ 22 ^^^r{R, V, u)cos{ye) \ e~^'^*dw, (20) 

^ J-oo U=o J 

where is the Neumann factor equal to 1 if = and to 2 otherwise. 

The series in (20) is non-absolutely convergent [10]. Its convergence in u depends on u 
and hence so does the truncation order which is supposed to satisfy the criterion 

\Ur{R, Z^m,^)| 

where e is the error bound. It is, however, known that usually the sum in question is 
mainly due to the first terms of the sum. For the frequency band used in our case, it was 
verified that taking a fixed truncation order Vjn such that keeps about 200 terms of the 
series ensures that the committed error does not exceed 10~^. 



(21) 



4-2. Results of numerical simulations 

Numerical simulations have been implemented for the opaque aluminium rod and the 

optically absorptive NG5 coloured glass rod, both of the diameter 5 mm. The values of 
their elastic and thermal constants are listed in Table 1 (the notations are explained below 
(1)). Assuming the laser wavelength 355 nm, the absorption length for NG5 coloured glass 
is = 0.8 mm. The typical values of the laser beam-width b = 0.15 mm and of the pulse 
duration 5 ns are taken. The radial component of the surface displacement Ur{R,9,t) is 
supposed to be detected at the point (r, 6) = (R, n) which is diametrically opposite to the 
generation point {R, 0) . Fourier transform in u> has been performed with the numerical 
parameter 5 = 0.035, sec (20). 

The results of numerical simulations are presented in Figs.l and 2. The wave arrivals are 
marked by the dashed lines. They are identified and denoted according to the nomenclature 
proposed in [4] on the basis of a ray trajectory analysis. The notation uL and uT 
corresponds to the arrivals of longitudinal and transverse bulk modes, travelling along the 
straight diameter path and undergoing n normal reflections from the cylinder surface at the 



Table 1. Physical properties of Al and NG5 glass used in simulations. 





P [g cm ^] 


A [GPa] 


^l [GPa] 




Cp [J kg-i K-i] 


Al 


2.69 


56 


26.5 


2.5e-6 


902 


XG5 


2.31 


17.4 


22.9 


G.50-6 


700 



edge points. The notation nP and nS also corresponds to, respectively, the longitudinal 
and transverse bulk modes, but these ones arrive at the detection point after a broken path 
resulting from n oblique reflections without modal conversion. The wave arrivals, involving 
modal conversion at oblique reflections, are denoted by mPnSr, where the subscript r is 
the number of round trips about the center made within the given path. The Rayleigh- 
wave arrival is denoted by Rr- The notation Hr is used for the wave arrival due to the 
head wave, i.e., due to the transverse mode radi ated by the s kimming longitudinal wave in 
the direction of the critical angle (pc = arcsin ^/2fi/ (A + 2/x). 

Fig.l compares two simulations of the laser-generated acoustic response Ur{R,'TT,t) for 
the opaque aluminium cylinder. The upper curve is computed by means of the developed 
procedure reduced to the model case such that assumes a radially concentrated buried 
source tending to the cylinder surface or, equally, a vanishingly small optical penetration 
length (see the closing remark in § 3.3). The lower curve is computed by means of Scruby's 
model of the surface dipole source [5] . Precise agreement of the arrival times and waveform 
shapes of both signals is observed. 

Fig. 2 presents the numerical results for the optically absorptive NG5 coloured glass 
cylinder. The upper curve shows the acoustic response Ur{R,7r,t) computed from Eqs. 
(18)-(20) for the radially distributed source, which takes into account marked enough 
optical penetration of the given material {aR 3). This curve is again compared with 
the 'would-be' signal calculated for this case on the basis of Scruby's surface-source model 
(lower curve) . It is seen that while the arrival times for both signals are obviously in good 
accordance, there is an evident discrepancy in the shape of the waveforms. For instance, the 
first longitudinal- wave arrival L on the upper curve has a characteristic bipolar broad shape, 
which is not observed on the lower curve. Another striking dissimilarity is that the bulk- 
wave and surface-wave arrivals on the upper curve have more or less similar amplitudes, 
which is in contrast to the predominant peaks of the first Rayleigh-wave arrivals R and Ri 
on the lower curve. Altogether, the difference between these two curves highlights the effect 
of optical penetration, which Scruby's surface-source model is not intended to describe. 

5. Conclusions 

Acoustic response of an isotropic cylinder to a laser line pulse has been simulated in 
the framework of 2D semi-analytical model, which takes into account optical penetration 
into the bulk of optically absorptive cylinder material. The waveforms calculated for the 
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Figure 1. Acoustic response of the aluminium rod computed for the model of a buried 
source tending to the surface (upper curve) and for Scruby's model of a surface source 
(lower curve). The notations for wave arrivals are explained in text. 



aluminium rod with a vanishingly small penetration coincide, as expected, with those 
obtained from Scruby's model of a surface dipole source. At the same time, calculation for 
the glass rod with large enough penetration length has demonstrated an essential difference 
of the waveform shapes relatively to the signal calculated for the same material by way of 
the surface-source model. 

The semi-analytical procedure developed in this paper paves the way to quantitative 
analysis of experimental results of laser ultrasonics generation in cylindrical samples made 
of optically absorptive materials. First experimental data have recently been obtained and 
will be reported elsewhere. 
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Figure 2. Acoustic response of the optically absorptive NG5 coloured glass rod. The 
upper curve is calculated for the radially distributed source taking into account optical 
penetration. The lower curve is calculated from Scruby's model of a surface source for 
opaque materials. 
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